Finite size Berezinski-Kosterlitz-Thouless transition at grain boundaries in solid 4 He 

and role of 3 He impurities 
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We analyze the complex phenomenology of the Non-Classical Rotational Inertia (NCRI) observed 
at low temperature in solid 4 He within the context of a two dimensional Berezinski-Kosterlitz- 
Thouless transition in a premelted 4 He film at the grain boundaries. We show that both the 
temperature and 3 He doping dependence of the NCRI fraction (NCRIF) can be ascribed to finite 
size effects induced by the finite grain size. We give an estimate of the average size of the grains 
which we argue to be limited by the isotopic 3 He impurities and we provide a simple power-law 
relation between the NCRIF and the 3 He concentration. 



The report of a Non Classical Rotational Inertia 
(NCRI) in solid 4 He 0, 0] has opened an intense de- 
bate in the physics community about its possible "super- 
solids" (SS) nature. Although the observation of NCRI 
has been confirmed by other groups [!, 0, 0], its phe- 
nomenology presents strong discrepancies with a simple 
supersolid phase, so that the precise origin of this phe- 
nomenon is still unclear. On one hand the supposed SS 
transition appears to be anomalously broad in tempera- 
ture [ESS [ill. On the other hand the NCRI strongly 
depends on the external conditions [1, 0, 0, 0, 0] • I n Refs. 
[HHl for instance annealing was shown to reduce and even 
to make disappear the NCRI fraction (NCRIF). More- 
over, using a different set-up, in Ref. [P| it was shown 
that the mass flow, associated with a SS phase, occurred 
only in the presence of grain boundaries (GB), and it 
was absent when GBs were not detected. This observa- 
tion gives rise then to an alternative hypothesis to the SS 
phase, namely, that a liquid phase is confined at the GBs 
and that the mass flow is related to superfluidity of the 
liquid component, similarly to a Rollin film. Partial wet- 
ting of GBs was experimentally observed in Ref. [lj| , and 
the possible superfluid (SF) ordering was investigated in 
Refs. Ill, LL2|, |l3( . Interestingly enough, a change of the 



shear modulus has also been observed at the low tem- 
perature, with a similar dependence on annealing and on 
3 He concentration as the NCRIF [3]. The connection 
between these two quantities is thus worth to be further 
investigated. 

Quite puzzling is also the dependence of the NCRI phe- 
nomenology on the 3 He concentration X3. The first re- 
port of NCRI [H was observed in commercial 4 He, which 
contains generally a low concentration X3 ~ 0.3 ppm. 
Further investigations showed that the critical tempera- 
ture increases monotonically with X3 whereas NCRIF in- 
creases with X3 only up to an optimal doping at X3 ~ 300 
ppb, after which the magnitude decreases [15| . 

In this paper we propose that, due to the strong con- 
finement on the grain boundaries, the SF transition of the 



premelted liquid component can be described in terms of 
a two dimensional Berezinski-Kosterlitz-Thouless (BKT) 
superfluid transition where the grain size gives rise to fi- 
nite size effects. We propose also a simple model where 
the concentration of 3 He impurities rules the grain size, 
and hence the finite size effects. We show that this frame- 
work can explain in a natural way, for X3 < 300 ppb, the 
broadness of the SF transition and the dependence of the 
NCRIF on the 3 He impurity concentration. 

In the following we shall model the polycrystal 4 He 
samples in terms of spherical grains with radius R and 
probability distribution function P(R). Premelting ef- 
fects, as discussed in Refs. 11, 12, 13, lij], is expected 
to give rise to a thin liquid film with thickness d. Partial 
wetting T3, 12 1, reducing the liquid amount covering the 
grain, can be also considered but it will not change our 
results. We assume that d <C R, so that the liquid he- 
lium system confined on the GB surface can be regarded 
as two-dimensional. This kind of model was employed by 
Kotsubo and Williams (KW) to explain the behavior of 



SF 4 He films on different substrates [17|. The important 
ingredient within this context is that the size R of the 
grain provides an intrinsic finite size cut-off which makes 
the BKT transition to be smooth. The broadness of the 
NCRI transition can be thus employed to estimate the 
average size i?o of the GBs. 

The BKT self-consistent equations on a spherical ge- 
ometry were discussed in Ref. [Tt} by KW. We can define 
the energy Uo(9) of an isolated vortex-antivortex pair at 
angular distance 8, in units of k^T, as 



U (9) = 2U C 
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where A'o = h 2 a® / 'm? k-gT is related to the bare areal SF 
density a®, m is the 4 He mass, U c is the vortex core en- 
ergy, and C = ciq/R is the minimum vortex-antivortex 
angular distance given by the vortex core size a^. In the 
presence of screening effects due to vortex pair polariza- 
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tion, we can generalize Eq. |T]) as 

where the static dielectric constant eo(6) can be evaluated 
in a self-consistent way as: 



1 



47T 3 /V 



d9'6' 2 sin & exp[-U(6')}. (3) 



Eqs. (HI)-© can be evaluated self-consistently for all 
9 < 7r to obtain the observable SF density a s (T) = 
Ug/e(8 = 7r) as function of temperature. Eqs. {I}-© can 
be also generalized to the dynamical case by introducing 
the dynamic dielectric constant e(9,u) which depends on 
the diffusion constant of the vortices D and on the an- 
gular frequency uj through the parameter rjj = \J2D /uj 



In the physical range of the experimental 
setup r_o/ao ^> 1 and the evaluation of the SF density a s 
in the dynamical regime is practically indistinguishable 
from the static one. The introduction of the dynamical 
analysis permits however the evaluation as well of the 

" mm. 



change of quality factor A[Q~ 

Let us now apply the above analysis to our polycrys- 
tal spherical-grain model. For a single grain of size i?, 
assuming d <C R, we can estimate the temperature de- 
pendent NCRIF n s {T,R) as 



n s (T,R) 



4irR 2 a s (T) 
4npR 3 /3 



3cr s 
~p~R 



(4) 



(where p is the solid 4 He density), and a mean n s (T) = 
J dRP(R)n s (T, R). In the following we shall show that 
the SF temperature profile is mainly ruled by the mean 
grain size value Rq- In this case we can roughly esti- 
mate the zero temperature NCRIF n s ~ Sa^/pRo. Note 
also that in the Ro ^ an limit the areal SF density 
a® is roughly proportional to the Berezinski-Kostcrlitz- 
Thoulcss temperature Tbkt where cr s (T) drops to zero. 
Within this context thus the effect of finite GB size R 
is mainly to give rise to a significant broadness of the 
transition whereas the temperature position of the drop 
is only weakly affected. If we assume a® ~ c s (Tbkt) = 
(2m 2 /7r/i 2 ) fee Tbkt, and we estimate Tbkt from the tem- 
perature T50 at which n s (T) drops to its 50 % value of 
n s {T = 0), we can thus obtain a free fitting parameter 
estimate of the GB size: 



Ro 



6m 2 k B T B KT 
irph 2 n s 



(5) 



We would like to stress that, because of the simplicity 
of this model and of the slight approximations in the 
estimates of Tbkt and Eq. ([5]) is simply meant to 
give the order of magnitude of R . 

We now apply our model to the specific case of NCR! 
in solid 4 He. We consider a 4 He density p = 0.2 g/cm 3 , 
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FIG. 1: (Color online) Comparison between experimental 
data for the 133 ppb TOP sample (symbols) [l5| and our the- 
oretical analysis. Also shown is the BKT curve for R — 00. 



which corresponds to a molar volume ~ 20 cm 3 and to a 
pressure 41 bar. We set also typical values for the vortex 
core size a — 50 A and energy U c = 2.5K [2(j. We 
consider for the moment two extreme probability distri- 
bution functions, namely a single value P(R) = S(R—Rq) 
and a flat P(R) = 1/R Q for R - AR <R<R + ARq, 
with ARo = Rq/2. With this choice of parameters, the 
overall profile of the NCRIF n s (T) is uniquely deter- 
mined by the only two free parameters, namely i?n and 
rj°, where cr°/Ro rules the magnitude of n s at T = 
while Ro is related to the broadness of the SF transition. 

In order to show the feasibility of this approach to 
reproduce the experimental results, we compare in Fig. 
[T] the n s vs. T data for the 133 ppb 3 He Penn State 
(TOP) sample [l5| with our best fit, which gives a® = 
0.26 x 10~ 9 g/cm 2 and Rq = 130 A, and with a infinite 
size BKT transition for large grains R 3> an. Note that 
there is only a slight difference between the (5-like and 
the flat P(R). As a matter of fact, we have checked 
that different distribution functions P(R) do not affect 
qualitatively our results, so that from now on we shall 
consider for simplicity a simple 5-like distribution. The 
nice agreement between the experimental data and our 
results suggests that the NCRI broad transition is not 
related to inhomogeneities of the samples but it stems 
from finite size effects due to the finiteness of the He 
grain. It is important to underline that, while the fitting 
procedure gives a refinement of the GB size, the order 
of magnitude of R is essentially given by the relation 
©, so that the good agreement in the broadness of the 
transition is not a result of the fit but it can be considered 
as an independent check of the validity of our analysis. 

The estimate size Rq — 130 A of the grains can appear 
quite puzzling especially considering that 4 He at these 
conditions is thought to solidify in a polycrystal form 
with macroscopic size of grains. However, on one hand 
the crystallographic evidence cannot esclude mosaics of 
small-angle grain boundaries. On the other hand, the 
presence of such small grains can naturally account for 
the high sensitivity to annealing and to preparation and 
freezing procedures, and the similarities with the NCRI 
phenomenology in Vycor glasses [2l|. Nevertheless, the 
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FIG. 2: (Color online) Comparison between theory and ex- 
perimental NCRIF and A[Q _1 ] data taken from Ref. fl5l ] for 
the 1 ppb TOF (upper) and CP (lower) samples. 

physical origin of such an extremely small grain size and 
its dependence on 3 He concentration requires to be ex- 
plained. We address this issue in the second part of this 
Letter, where we relate the origin of GB to the presence 
of liquid melted bubbles induced by 3 He impurities. 

In order to gain a further insight on this point, let us 
discuss before the experimental dependence of the NCRI 
phenomenology on the He 3 concentration and on the 
growth/measurement condition. In particular, we ana- 
lyze two samples [l5[ in the highly dilute limit £3 ~ 1 
ppb, one belonging to the University of Florida (TOF) 
family, grown with the blocked capillary method, and one 
still belonging to the Penn State University but grown at 
constant pressure (CP). For both these samples measure- 
ments of the change of the quality factor A[Q _1 ] were 
also available [lit - 

The first thing to observe is that, despite both sam- 
ples are at same £3, the signals are significantly different 
(Fig. [2]). For the TOF sample, the value of n s is al- 
most one order of magnitude smaller than for the CP 
sample. Its variation as a function of the temperature 
is also different. In fact, while the n s presents a quite 
smooth behavior in the CP sample, the drop of n s at 
T ~ 0.03 K is much sharper followed by an additional 
tail. Similar features are observed in the behaviors of 
the quality factor A[Q _1 ] which is quite broad in the CP 
sample compared to the sharp peak at T ~ 0.03 K in the 
TOF sample (note that the longer tail in the n s is not 
observed in A[Q _1 ], suggesting that this features in the 
n s (T) behavior could have a spurious origin). 

We can see now that our model permits to understand 
in a very natural way also the X3 dependence of the NCRI 
phenomenology. We first note that, although the two 
samples have different n s (T = 0) of a order of magnitude, 
their Tbkt, defined for instance as Tbkt = T50, are of 
a similar magnitude. A simple analysis, using again Eq. 
([5|), would point out thus an average radius Rq of the 
GBs in the TOF sample one order of magnitude larger 
than in the CP one. As a consequence, finite size effects 



are expected to give rise to a much sharper transition 
in the TOF than in the CP case, in agreement with the 
experimental observation. These simple considerations 
are corroborated by a more detailed analysis. In Fig. [2] 
we show our best fits for both the TOF and CP samples, 
compare with the experimental data. Estimates for Rq 
and (To, in these cases, are respectively Ro = 790 A, a® = 
0.105 x 10~ 9 g/cm 2 , and R a = 100 A, ct° = 0.216 x 
10~ 9 g/cm 2 . For the evaluation of the quality factor 
A[Q _1 ] we have used respectively rjj = 1.2 x 10 4 A and 
rjj = 0.3 x 10 4 A, which give rjj/R 1. We remind 
that in such quasi-static regime the parameter r^> simply 
acts as a scale factor on A[Q _1 ] while it does not affect 
significantly the NCRIF. The reliability of such limit is 
confirmed by Ref. [f| where the magnitude of A[Q _1 ] 
was shown to be sensibly affected by changing the set-up 
frequency while the NCRIF was essentially untouched. 

We can see that the finite GB size BKT theory can nat- 
urally account not only for the broad NCRIF transition 
but also for the experimental height and broadness of the 
A[Q -1 ] factor. This latter point is not surprising since in 
the finite size BKT framework, the broadness of the n s 
transition and of the A[Q _1 ] peak are related. Although 
our theory can account for the difference of the signals 
in a simple manner, the long tail of n s for T > Tbkt 
in the TOF sample remains unexplained. As mentioned 
above, however, it should be noted that this feature of the 
NCRIF has no counterpart in A[Q _1 ], which is sharply 
located at T Tbkt UM- This suggests that the long 
tail in n s (T) might be due in this sample to different 
physics not related to the supcrfiuid GB transition. 

After assessed the robustness of our analysis in ex- 
plaining the NCRI phenomenology as function of X3 and 
growth conditions, in the last part of this Letter we ad- 
dress the origin itself of such a^-dcpcndcncc. In particu- 
lar we propose that the main effect of 3 He impurities is to 
provide an intrinsic maximum length scale for the growth 
of grains, favoring thus the presence of grain boundaries, 
and hence to sustain a supcrfiuid NCRI on the GB sur- 
faces. It should be clear, on the other hand, that ad- 
ditional sources of disorder, favoring the GB formation, 
can be present, as shown by the different phenomenol- 
ogy for differentgrowth conditions and by the annealing 
dependence [H, @|. We assume for the moment that the 
only source of GB is the presence of 3 He impurities. This 
assumption is probably valid for the TOF samples, which 
show the smallest NCRIFs for similar X3, and it is cor- 
roborated by the sharp drops of n s at Tbkt, suggesting 
quite large grains. We consider, for dilute 3 He concen- 
trations 2:3 < 1 ppm, a uniform distribution of the 3 He 
impurities within the sample, with an average distance 
d3Ho between 3 He atoms dsjjc ~ £3 1//3 . In Ref. (ljj it was 
shown that 3 He impurities, due to the stronger quantum 
fluctuations of zero point motion, induce a local melting 
of the host 4 He even at low temperature much smaller 
than the bulk 4 He melting. The presence of local liq- 
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FIG. 3: (Color online) (a) schematic sketch of GB structure 
of 4 He samples in the presence of 3 He impurities. Light blue 
areas represent solid 4 He, white spots 3 He impurities and dark 
blue regions premelted 4 He at the GBs and around 3 He impu- 
rities, (b) n„ vs. X3 plot of different samples (symbols) taken 
from Ref. [l5| compared with our model prediction [Eq. ] 
(solid line), which is expected to fail in the gray region where 
R > 2 — 3ao and the mean-field theory breaks down. 



uid spots, which we remind survive also at virtually zero 
temperature, around the 3 He impurities poses a strong 
constraint on the growth of the grain size in the freez- 
ing process. The average distance c?3 He provides thus the 
maximum length scale for the growth of grain size and, 
in the last analysis, an average value of the grain diame- 
ter 2i? — He (see Fig. for a sketched picture). This 
simple estimate gives for instance Rq ~ 1600 A for X3 = 1 
ppb and i?o — 313 A for £3 = 133 ppb, 2-3 times larger 
than the actual fit estimates. Taking into account that 
these figures are purely indicative since other sources of 
disorder are always present further limiting the size of 
the grains, such estimates are not bad and provide the 
order of magnitude of the grain size. The most convinc- 
ing probe of such picture is the power-law dependence of 
the behavior of n s as function of X3, n s « x^ 3 . Assum- 
ing that a p s fraction of the thin liquid film of thickness d 
at the GBs undergoes a superfluid BKT transition, and 
approximating the liquid 4 Hc density with the solid one, 
we have an areal SF density at the GBs a s = p s dp, and, 
from Eq. Q , 



6p s d/d 3lie 



6p s dxl /3 /a, 



(6) 



where we have assumed d <^ Rq and i?o ~ d,3n c /2. where 
a = 3.2 A is the solid He-He distance at the pressure 
here considered. Such power-law behavior, assuming for 
instance d = 1.2 A and p s = 0.01, is shown in Fig. [Hb 
in qualitative agreement with the experimental data for 



the TOF samples and the lowest X3 TOP samples [15 1 
where other sources of disorder are thought to be small. 
Note that the unknown quantities p s , d, in the log- log 
plot of Fig. [3Jd, determine only the vertical off-shift of 
the logn s — logX3 behavior but not its slope which is 
uniquely determined by geometrical considerations. In 
CP samples on the other hand additional limiting mech- 
anisms on the grain size are probably operative conceal- 



ing the 3 He effects. This is in agreement indeed with the 
larger values of n s and with the broader drops of NCRIF 
@, Ull- The validity of this analysis is in addition lim- 
ited by the mean-field character of our approach. For 
x 3 ^ 10 3 ppm the average size of the grains is predicted 
to be Rq < 150 A, only 2-3 times larger than the vor- 
tex size do = 50 A. In this regime the reliability of the 
mean-field analysis is strongly questionable. 

In conclusion, in this Letter we showed that the com- 
plex NCRI phenomenology is fully compatible with a 
superfluid Bcrczinski-Kostcrlitz-Thouless transition in- 
duced in the thin premelted liquid film at the grain 
boundaries. Note by the way that the two dimensional 
BKT transition is characterized by the lack of specific 
heat anomalies, in agreement with the heat-capacity 
measurements in solid 4 He [22], [23[ [24[ ■ Both the temper- 
ature and 23 dependence are shown to be ascribable to 
finite grain size effects. We propose also a simple picture 
where 3 He impurities are directly related to the maxi- 
mum size of the grains and we predict a simple scaling 
relation between the NCRIF n s and the 3 He impurity 
concentration. 

We acknowledge useful discussions with J. Beamish, H. 
Alles, A.V. Balatsky. M. Holzmann, S. Balibar and F. 
Caupin. 
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